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Abstract 

Using the recently released Union2 compilation with 557 Type la su- 
pernovae, the shift parameter of cosmic microwave background given 
by the WMAP7 observations, and the baryon acoustic oscillation mea- 
surement from the Sloan Digital Sky Survey, we perform the anal- 
ysis on the 1-loop Yang-Mills condensate (YMC) dark energy model. 
The analysis has been made for both non-coupling and coupling mod- 
els with flrnO and wq being treated as free parameters. It is found 
that, xizn = 542.870 at QmO = 0.2701 and wq = -0.9945 for non- 
coupling model, and Xmin — 542.790 at 7 = -0.015, Q^o = 0.2715 and 
wq = -0.9969 for coupling model. Comparing with the ACDM model, 
the YMC model has a smaller xLin^ but it has greater values of the 
Bayesian and Akaike information criteria. Overall, YMC is as robust 
as ACDM. 
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1 Introduction 



The discovery of the current accelerating expansion of the Universe via the 
observations of the SNIa |I] with follow-ups [2], further supported by the ob- 
servations on CMB [21, S] and on large scale structure (LSS) fSj, has brought 
a most challenging problem on both cosmology and physics. An abundance of 
literature has been devoted to the issue. To interpret the cosmic acceleration, 
one may either go beyond the general relativity (GR) and seek an alternative, 
such as the effective gravity [B] , or stay within the framework of GR and simply 
attribute the acceleration to some mysterious dark energy (DE) as the driving 
source. A positive cosmological constant A is the simplest candidate of DE, but 
it is plagued with the coincidence problem [7]. Various dynamic models have 
been proposed address this issue p]. One class of models is based upon some 
scalar field, such as quintessence [9], k-essence pO], phantom |IT], quintom fi2\ , 
etc. In our previous works 1131 [HI [El [IS] , we have developed a vector field type 
of dynamic DE model, in which the renormalization-group improved effective 
Yang-Mills Condensate (YMC) serves as the dynamical DE. This has bee 
highly motivated by the great success of Yang-Mills fields as a corner-stone 
of particle physics, which mediate interactions between fundamental particles 
and define the vacuum structure. It has been demonstrated that YMC DE 
model has the following properties desired for a dynamical model: being able 
to solve the coincidence problem naturally, giving an equation of state (EoS) 
w that cross —1 smoothly in the coupling case, having the dynamic stability, 
and alleviating the high-redshift cosmic age problem. Moreover, these proper- 
ties are retained by YMC models even with the increase of order of quantum 
corrections, i.e, all 1-loop pL4j, 2-loop [T5], and 3-loop models have ex- 
hibited the same dynamical behavior. For a detailed description of the 1-loop 
YMC dark energy model, including its theoretical basis and the dynamics, see 
Ref.[14]. Other models using spin-1 fields as a candidate for the dark energy 
can be found in Refs.|I7]. 

Any DE model has to be confronted with observations. Recently, there are 
some significant updates in the observations of SNIa and CMB. The Union 
SNe la compilation [[TD] enlarged by the CfA3 sample [20] has been recently 
updated, including a few SNe la of high redshifts, and contains 557 SNe la, 
forming the Union2 compilation [21]. Not only the number of SNe is sub- 
stantially increased, but also the range of redsifts is extended, crucial for de- 
termining the evolution of dynamical dark energy. Besides, the 7-year data 
of Wilkinson Microwave Anisotropy Probe (WMAP) has given an improved 
determination of cosmological parameters [22], in combination with the latest 



distance measurements [23] of the Baryon Acoustic Oscillations (BAO) ||24] in 
the distribution of galaxies of SDDS and 2dFGRS. Motivated by these, it is 
natural for us to expand our work ||TB] to further constrain the YMC model 
with these observational datasets. In this work, we shall carry out a statistical 
analysis, utilizing the combination of SNIa, CMB, and LSS data, to constrain 
1-loop YMC model for both non-coupling and coupling cases. The examina- 
tion is much more refined than the previous work |[TB] in performing statistics. 
Furthermore, employing the Bayesian information criterion(BIC) [27] as well 
as the Akaike Information criterion(AIC) [28], we shall also perform extra sta- 
tistical examinations for models with different number of parameters. The 
resulting statistics shows that, comparing with the ACDM model, the YMC 
model is still a robust dynamical dark energy model. The unit system with 
c = h = 1 is used in this paper. 



2 The 1-loop YMC model 

We consider a spatially flat {k = 0) Robertson- Walker (RW) universe, whose 
expansion is determined by the Friedmann equations 

H^ = ^iPy + Pm), (1) 
d 4:7tG . , , 

- = :^{py + ^Py + Pm), (2) 

where H = ^, py and pm are the energy density of the YMC and the matter 
(including both baryons and dark matter), respectively, Py is the pressure of 
the YMC. For simplicity, the radiation component is neglected since its con- 
tribution is very small in the matter dominated era under consideration. The 
energy density py and the pressure py of the 1-loop YMC are given by [221 



Py = lb,,\y + l)ey, (3) 

Py = Ib^ily - (4) 

where k, is the renormalization scale of dimension of squared mass, b = -^^^ 
for the gauge group SU{2) without fermions, y = ln \ E'^/k,'^\ with being the 
squared electric field of YM condensate . When one requires that py be the 
dynamical dark energy and its value at z = be equal to ~ 0.73pc, one finds 
K^/^ ~ 5 X 10-^ eV [II]. The equation of state (EoS) for the YMC is 

w = — = -. (5) 

Py 32/ + 3 ^ ^ 



The dynamical evolutions of the YMC dark energy and the matter are given 
by 

Py + ^-{Py+Py) = -^Py. (6) 

(Jj 

Pm + 'i-pm = Tpy, (7) 

Gj 

where F is a model parameter, representing phenomenologically the coupling 
between the YMC and the matter components. For F > 0, the interaction term 
Tpy is the transfer rate of the YMC energy into matter, whereas F < means 
that the matter transfers energy into the YMC. For computing convenience, 
Eqs.([Ij), (ED and ([Tj) can be recast into: 

= ^JI^h\ (8) 

dx {l + y)ey 
dN = 

^ = - ^ + ^ (-in) 

dN ^{2 + y)h 2 + ^ ^ 

where x = pm/\bti^, N = Ina(t), h = + y)ey + x, and 7 = F/ ^/47TGlmP^. 

The set of dynamic equations is completely determined by the initial values 
Xi and yi for the noncoupling case, plus 7 in the coupling case. In actual 
computation, the initial condition can be taken at a redshift Zi = 1092. For 
a given 7, each set of (x^, yi) yields a corresponding set (xq, yo) at z = 
as the outcome from solving the dynamic equations, and hence EoS wq = 
{yo — 3)/{3yo + 3) and the matter fraction fl^Q = xo/{xq + {yo + l)e^°). This 
point of our treatment differs from those models like XCDM [30], where wx or 
QmO in some model, was put in by hand as parameters, instead of following from 
dynamic equations. It should also be mentioned that we will not do statistical 
examination with the Hubble parameter Ho = 100 h km see-^Mpc-' in this 
paper, (xq, yo) can be taken as two independently adjusted parameters for the 
statistical examination. 



3 Constraints from SNIa, BAO and CMB 

Below, we confront the YMC model with the latest observational distance 
modulus Pobs{zi) data of 557 SNIa [21], the BAO measurement from the Sloan 
Digital Sky Survey (SDSS) [23] and the shift parameter of CMB updated by the 



7-year WMAP observations [22] • The theoretical distance modulus is defined 
as 

lith{z) = 5\og^QDL{z) ^ llQ (11) 

where 

/io = 42.38- 51ogio/i, (12) 

and 

is the luminosity distance in a spatially fiat Universe, actually independent of 
the Hubble constant Hq. For our model, 



^(z) = Wq.„o^ + (1-0„.o)^, (14) 

V Pm{0) Py{0) 

which depends upon the adjusted parameters (ilmO , 7). Here E{z) is 
an implicit function for YMC through Pm{z) and Py{z) as the solution from 
the dynamical equations (ISjl-ffTUj). In doing the statistical examination in the 
following, an ensemble of the solutions of the equations are generated, each of 
which correspond to a point in the grid of (n^o ? 7). This requires much 
more computing time than the models like XCDM with an expression for E[z) 
containing parameters explicitly [301 EI] • 

Note that, to reveal the possible dynamical property of dark energy, we 
choose the EoS wq as free parameter, which is of dynamical nature and re- 
fiects the second order time derivative a{t). This is pertinent since the Union2 
dataset provides SN la with higher redshifts, better for constraining the evo- 
lutional property of a DE model. For examining a dynamical DE model, this 
has the advantage to the choice of the Hubble parameter h in Ref.[IB]. More 
importantly, the constraints on cosmological parameters by the BAO measure- 
ment from SDSS and 2dF Galaxy Redshift Survey (2dFGRS) has an intrinsic 
degeneracy with (^1^0, h). Therefore, one can not arrive at reliable estimate for 
both and h simultaneously [23, 121] • As a matter of fact, if one would want 
to determine h with sufficient accuracy, one has to go beyond and to employ 
some extra data, such as the HST's key project [25] and SDSS for the history 
of the Hubble parameter [2S], etc. 

For the SNIa data, the corresponding estimator is constructed as: 

557 2 

xL(p;w) = E '''°'°^''^"2'"'^'''' ■ (15) 

— (77 



where p stands for a set of parameters, such as {^mo^ wq, 7), fiobs{zi) is the 
observed value of distance modulus for the ith supernova, and cr^ is the cor- 
responding la error that can be found via Ref.|[2I]. The nuisance parameter 
fiQ can be analytically marginalized over [31], so that one actually minimizes 
Xsn{p) instead of X^sn{P'^ f^o)- The minimization with respect to fiQ can be 
made simply by expanding the of Eq. ( TTHI ) with respect to jj.o as 

xIAp) = Mp) - 2/ioi^(p) + f^lC, (16) 



where 



A{p) = ~ ^^th{z^]^o = 0,p)]^ ^^^^ 



=1 
557 



[l^ohs{Zi) - llth{Zi] 1^0 = 0, p)] 
(7. 



1=1 ' 



C = Ei- (19) 

1 = 1 « 

Evidently, Eq.([16j), as well as (fT5|). has a minimum for fiQ = B/C at 



Xhip) = A{p) - (20) 

Since XsN^mtn = XsN^Tmn^ lustcad of minimizing Xsn one can minimize xliv 
which is now independent of the nuisance parameter /ig. Notice that this ana- 
lytical marginalization over fiQ means that the Hubble parameter h is effectively 
marginalized over. 

Next, the distance parameter A of the measurement of the BAO peak in the 
distribution of SDSS luminous red galaxies is defined as |[2i 



1 dz 



2/3 



(21) 



.Zb Jo E[z) 

with Zh = 0.35, the redshift at which the acoustic scale has been measured. 
The observations [2^ give 

Aobs = 0.469(n,/0.98)"°-^^ ± 0.017, (22) 



where Ug is the primordial spectral index and WMAP7 data [22] yields the 
updated value = 0.963, while it was = 0.960 by WMAP5 [4]. The 
corresponding of BAO is given by 

xIao = (23) 
"a 



where a a = 0.017. 

Finally, the shift parameter R from CMB is defined as 

with Zrec being the redshift of recombination. WMAP7 [22] gives Zrec = 
1091.3 ± 0.91 and Robs = 1-725 ± 0.018, while it was Zrec = 1090.0 ± 0.93 
and Robs = 1-710 ± 0.019 by WMAP5 H]. The of CMB is 

XCMB = 2 ' 

(25) 

with aji = 0.018. As usual, assuming these three sets data of observations are 
mutual independent, and the measurement errors for each set are Gaussian 
with the likelihood function of the form 

jC oc e~^'/l (26) 

The three sets data are combined by multiplying the likelihoods, and the com- 
bined is given by 

= xIn + x\ao + xImb- (27) 

Variations of values of the parameters (ilmO, wq^ 7) yield respective values of 
X^. A minimal has been found, corresponding to a maximal which would 
be favored by the observations. 



4 Results 



For the non-coupling 7 = case, there are two parameters Q^o and wq. In 
our present computation of the confidence level of ^,720, the procedure goes as: 
First we compute x^(^?7i05 wq) and the Likelihood C = C (r^mO, wq) as functions 
of both n,„o and wq, from which, by a standard searching procedure, follows 
the resulting minimum Xmin = 542.870 at Q^o = 0.2701 and wq = —0.9945. 
Next we integrate C (0^0, wq) over wq and derive the Likelihood C (Qmo) as a 
function of only Q,„o- The prior cut off wq = (—0.9999, —0.90115) is assumed 
here. Then we compute the confidence level of based upon the function C 
(Q^o)- The result is that = 0.2701+5j;[j;g at 68.3% CL(Confidence Level), 
Vt^o = 0.2701^i;:[jg5 at 95.4% CL. Since the non-couphng YMC model always 
has wq > —1 IH], it is found that the likelihood function C{wo) is mainly 
distributed in a narrow range close to = — 1, so that the computation of 
CL can only be done for an upper bound of wq. The result is that wq < 
-0.9773 at 68.3% CL, and wq < -0.9579 at 95.4% CL. The prior cut off 



^mo = (0.2001,0.3501) is assumed here. Figure 1 shows the details of the 
C.L. curves in the (r^mO, w^o) plane. We remark that, in obtaining the values of 
C.L. for a Q^no (or wq) at the maximum of £, the integration over wq (or flmO^ 
respectively) has been carried out. 

To make a comparison with the ACDM model as a sort of referring model, 
we have also done the computing and fitting of the ACDM model to the same 
observational data as above, yielding Xmin ~ 542.919 at Q^o = 0.2701. The 
confidence levels are the following: = 0-2701]^qo|35 at 68.3% CL, Q^o = 
0.2701 

-ao264 95.4% CL. We present the corresponding and hkelihood 
jC oc e^^ in figure 3. 

Based on the estimator only, one would draw a conclusion that the non- 
coupling YMC model is more favored in confronting this set of observations, 
as it has a slightly smaller xLin than the ACDM model. A similar result was 
obtained for non-coupling YMC model at 3-loop for an earlier dataset of the 
182 SNIa +CMB+BAO |I6] . This kind of better performance of YMC is due to 
its dynamical nature, say, its w{z) being a function of z, instead of a constant 
— 1 as in ACDM. In this regards, more observational data at high redshifts are 
desired to distinguish various DE models with different evolutionary behavior, 
by using certain diagnosis, such as the Statefinder and Om [32] • 




Figure 1: Non-coupling YMC: the C.L. curves in 0^0 — "'^'o plane based upon the joint study of SNIa, 
BAG and CMB. 

For the coupling model, there are three independently adjusted parameters 
7, flmO and wq. In order to search for possible DE models in an extended 
domain, in particular, we will allow the coupling 7 to take negative values. 
For each fixed value of 7, we minimize with Q^q and wq. We have done 
ten different values of 7, and the resulting values of Xmin corresponding 
sets of (7, flmO Wq) are listed in Table 1. It is seen that, as a function of 7, 
the minimum Xmin — 542.790 is attained at 7 = —0.015, Q^o = 0.2715 and 




Figure 2: Coupling YMC with 7 = —0.015: the C.L. curves in — wq plane based upon the joint 
of SNIa, BAO and CMB. 
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Figure 3: The and C of ACDM based upon the joint of SNIa, BAO and CMB 



wq = —0.9969. This Xmin slightly smaller than that of the non-coupling and 
of the ACDM models. Table 1 also tells that, although the Xmin depends on 7, 
the dependence is not very strong. 

Our previous work |[TB] did not search for the minimal with respect to 
the parameter 7 since it calculated only three different values of 7. 

Figure 2 shows the C.L. curves with 7 = —0.015 in the (QmO, wo) plane. 

Thus the YMC model with a negative coupling 7 < does a bit better than 
the latter two models in confronting the updated observations. We remark 
that this conclusion was partly hinted in the 3-loop YMC model, where the 
consideration was confined only to the limited region of 7 > |[16]. Note that 
7 < means a situation in which the matter is transferring energy into the 
YMC, and EoS wq will not cross over —1 IHj. Besides, in the far future when 
the scale factor is hundreds times the present one, the matter density pm will 
turn into negative, which would be a non-physical region. Thus, the model 
of a negative coupling is phenomenological and can not be infinitely extended 
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Table 1: For each given coupling 7, the values of Xmm corresponding best-fit parameters 

and Wo are listed. 



Model 


ACDM 


non-coupling 


coupling 


x'^ ■ 


542.919 


542.870 


542.790 


k 


1 


2 


3 




0.973 


0.975 


0.976 


ABIC 





6.277 


12.523 


AAIC 





1.951 


3.871 


Rank 


1 


2 


3 



Table 2: comparison of models 

into far future. In regards to the issue of coupling, we notice that a positive 
interaction is preferred in a study of dynamics of galaxy clusters with a coupling 
between dark energy and dark matter ||33] . 

The analysis is effective in searching for the best-fit values of parameters 
within a given model. But for models with different number of parameters, one 
would expect Xmm decreases as the number of free model parameters increase. 
For model comparisons in this case, one can use other kinds of criteria for 
model selection, such as BIC |27], AIC [2B], and Xmin/^of with the degree of 
freedom dof = N — whereas and k are the number of data points and 
the dimension (number of independently adjusted parameters) of the statistical 
model, respectively. The BIC is defined as BIC = —2 In jCmax + klnN , where 
J^max is the maximum likelihood. In the Gaussian case, Xmin ~ ~'^^^^max + 
constant, so that the difference in BIC is given by ABIC = Ax^in + ^• 
The AIC is defined as AIC = —2ln jCmax + 2A: , and the difference in AIC is 
given by AAIC = Ax^in + 2 A/c. We present the results for these three criteria 
in Table 2, which imply that the ACDM is still favored due to its simplicity, 
whereas the coupling model is ranked last as it has most parameters. 

In summary, using the x^ analysis alone in confronting to the updated obser- 



vational data of 557 SNIa +CMB+BAO, we find that the non-couphng YMC 
is favored over ACDM, and the coupling YMC with a negative 7 is the most 
favored. Taking into account of the dimension of a statistical model, ACDM is 
still simplest. Overall, YMC is as robust as ACDM. More observational data 
at higher redshifts in future are much desired to shed light on the issue of dark 
energy. 
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